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There is an urgent need to develop predictive indicators of the effect of species loss on ecosystem 
functioning. Body size is often considered as a good indicator because of its relationship to extinction risk 
and several functional traits. Here, we examined the predictive capacity of species body size in marine and 
freshwater multitrophic systems. We found a significant, but weak, effect of body size on functioning. The 
effect was much stronger when considering the effect of body size within trophic position levels. Compared 
to extinctions ordered by body size, random extinction sequences had lower multiple species loss effects on 
functioning. Our study is the first to show experimentally, in multitrophic systems, a more negative impact 
of ordered extinction sequences on ecosystem functioning than random losses. Our results suggest apparent 
ease in predicting species loss effect on functioning based on easily measured ecological traits that are body 
size and trophic position. 

The effects of biodiversity loss on the functioning of natural ecosystems are not easily predictable 1 , especially 
for complex food webs 2 . They might depend on the competitive dominance or the trophic position of species 
lost, interaction strengths between species, functional traits of both species lost and those remaining in the 
system, and relative control of both biotic and abiotic factors over ecosystem properties 3 . A major challenge is to 
derive predictive and easily measurable indicators of biodiversity effects on functioning accounting for their 
complexity. Recent studies suggested that the effect of single species loss could be predicted without complete 
investigations of system interactions 4 . 

Central to the metabolic theory of ecology, body size is linked to multiple biological rates such as growth, 
reproduction and mortality 5,6 . Measuring body size is an easy way to collapse co-varying traits over a single 
dimension, with no need to observe these traits directly 7 . In nature, species removal rarely happens randomly, 
contrary to simulations performed in past studies 8 . Larger species are often considered more vulnerable to 
extinction 9 than smaller ones because of low population densities, slow growth, high energetic requirements 
and overexploitation 10 . Simulations also showed that losses of larger species would have a greater effect on 
ecosystem functions, such as bioturbation, than smaller species 9 . 

We experimentally tested the predictability of body size for ecosystem functioning. We conducted experi- 
mental extinctions in multitrophic systems of up to 10 taxa in marine and freshwater mesocosms and assessed 
subsequent effects on ecosystem function. We predicted that 1) the effect of removing a taxon on the system will 
increase with its body size. Then, we compared the effect of random and non-random multi species sequential 
extinctions. We predicted that 2) non-random sequential extinctions (ordered by body size) will increase the rate 
of change in ecosystem functions compared to random extinction sequences. 

Results 

We found a significant positive effect of body size on marine multifunctioning (1st PCA axis only; adj. R 2 = 0.14, 
Fig. 1, Table 1). Body size had a negative effect on phytoplankton biomass (adj. R 2 = 0.30) in the freshwater 
ecosystem (Table 2). We also found a significant effect of body size on encrusted net production in the marine 
ecosystem once controlling for trophic position (adj. R 2 = 0.14). The addition of trophic position also increased 
the proportion of variance explained by body size in marine multifonctioning (adj. R 2 = 0.33; AAIC = 6.06) and 
phytoplankton biomass (adj. R 2 = 0.44; AAIC = 5.86). 

Asymmetrical ANO VA showed that the effect sizes of periphyton dry mass and marine multifunctioning were 
greater for extinction sequences ordered by body size, and the phytoplankton biomass lower, compared to 
random extinction sequences (Fig. 2). In fact, the effect of the first taxon removed (the intercept) was greater 
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Figure 1 | Effects of body size on ecosystem functioning. The proxies of ecosystem functioning illustrated here are (a) marine multifimctioning, (b) 
freshwater phytoplankton biomass and (c) marine encrusted algae net production. • indicate grazers and O indicate non-grazers. The y-axis refer to the 
effect size between values of functioning recorded in mesocoms with removals and values of functioning recorded in reference mesocoms (no removal). 



for ordered sequences in both systems, but then the effect of sub- 
sequent removal converged and the slope of the extinction-function- 
ing relationship was not different between the two sequences. We 
found no significant effect of the number of taxa removed on other 
functions. 

Discussion 

The experiment, performed with two different aquatic ecosystems, 
showed a stronger effect of large taxa extinctions on ecosystem func- 
tions (supportive of prediction # 1 ) . The removal of smaller taxa from 
the marine ecosystem, such as Littorina spp. and T. testudinalis bene- 
fited encrusted algae net production, while macroalgae net produc- 
tion increased when larger taxa such as S. droebachiensis were 
removed. Conversely, the removal of large taxa in the freshwater 
ecosystem decreased phytoplankton biomass. It is noteworthy that 
body size had predominantly affected functions directly related to 
consumption. This result might be explained by increasing strength 
of per capita trophic interactions with body size 4 . 

Although body size can be viewed as a measure aggregating several 
traits influencing ecosystem functioning 7,9,11 , it was not sufficient to 



explain alone the variability observed in functioning. The addition of 
the trophic position as a second explanatory variable revealed a 
stronger relationship between grazers body size and functioning than 
for the other taxa. The predictive capacity of the body size was there- 
fore much more convincing when combined to other traits related to 
trophic position. Extinction risk and its subsequent effects on func- 
tioning might therefore depend on the trophic position, with greater 
effects associated with the loss of large consumers 12,13 . 

We also found contrasting effects of random and non-random 
multi-taxa extinction sequences, in agreement with simulation stud- 
j es 9,i4-i6 j n me marme ecosystem, both sequences overall led to 
higher periphyton dry mass and multifunctioning, but non-random 
extinction sequences impacted more on the functioning. Conversely, 
the decrease in phytoplankton biomass following non-random 
extinctions occurred at lower values of biomass in the freshwater 
ecosystem. For both systems, the effect of the first taxon removal 
was stronger for the ordered than the random sequences, but then the 
slope of the extinction-functioning relationship was the same. 

As mentioned above, larger taxa have greater risks to get extinct 9,10 . 
Here, our results suggest that in addition, the consequences of losing 



Table 1 | Statistical results for the marine system, a) linear regression, b) ANCOVA with trophic position as predictor and c) asymmetrical 
two-way crossed ANOVA for the marine system with Number of taxa removed (no) and Sequence as main factors. The variability of the 
factor Sequence was then divided between 1 ) random versus non-random extinction sequence (Rvs NR) and 2) differences among random 
extinctions sequences (among random). P values in bold correspond to significant results 

Marine system 



Sources of variation 



Macroalgae net Encrusted algae net Multifunctioning Multifunctioning 
production production (I s1 axis) (2 nd axis) 



Periphyton dry mass 



(a) 


Adjusted R 2 


0.0432 




0.0000 




0.0985 




0.1383 




0.0000 






Body size 


2.1745 


0.1528 


0.5843 


0.45 1 8 


3.8396 


0.0613 


5.1726 


0.0318 


0.4449 


0.5109 


(b) 


Adjusted R 2 


0.0568 




0.3093 




0.1373 




0.3341 




0. 1 640 






Body size 


1 .4243 


0.2444 


2.8514 


0.1042 


5.1714 


0.0322 


9.9722 


0.0043 


0.0250 


0.8758 




Trophic role 


1.3587 


0.2552 


12.7837 


0.0015 


2.1259 


0.1578 


8.3510 


0.0081 


6.5560 


0.0172 


(c) 


Adjusted R 2 


0.2569 




0.4206 




0.4158 




0.5570 




0.4422 






No. taxa removed 


5.431 1 


0.0143 


8.6657 


0.0023 


0.631 1 


0.5434 


5.6244 


0.0127 


1.1036 


0.3531 




Sequence 
























Rvs NR 


8.7416 


0.0044 


0.1243 


0.7256 


1 . 1 647 


0.2848 


6.2287 


0.0153 


1 1.8585 


0.001 1 




Among random 


2.0768 


0.0522 


3.5819 


0.0019 


4.3657 


0.0003 


5.0189 


0.0001 


4.8956 


0.0001 




No. X Sequence 
























No. X (R vs NR) 


1.5384 


0.2231 


0.3216 


0.7262 


0. 1 004 


0.9046 


0.3542 


0.7032 


1.8618 


0.1643 




No. X (Among random) 


1.1490 


0.3343 


2.0241 


0.0257 


3.2733 


0.0004 


3.5915 


0.0002 


2.4692 


0.0060 
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these larger taxa might be more significant for the ecosystem func- 
tioning. The conjugation of these two factors indicates that we should 
be more concerned by the non-random character of species extinc- 
tions. To our knowledge, our study is the first to show experimentally 
the negative effects of ordered extinctions sequences in multitrophic 
systems. However, the loss or gain of a particular function, as the 
number of taxa removed increases, did not occur more rapidly fol- 
lowing a non-random extinction, contrary to our second hypothesis. 
It was previously shown that the complexity of interactions within 
ecological food webs makes biodiversity-ecosystem functioning rela- 
tionships idiosyncratic and almost unpredictable 2,17 . We indeed 
observed a much larger variability among random sequences. The 
analysis of identity effects in this system revealed that species inter- 
actions in food webs might balance each other and therefore prevent 
the observation of any general relationship between biodiversity and 
ecosystem functioning 17 . 

Our experiment reports the immediate effect of removing the 
largest organisms on ecosystem functioning. Our hypothesis was 
based on the idea that larger organisms have higher absolute energy 
consumption. However, metabolic theory also predicts that the effi- 
ciency of energy consumption should increase with species body 
size 6 . Consequently, over the long-term, we should expect that the 
loss of large species would be compensated by less efficient smaller 
ones. It was not possible to test for this second mechanism due to the 
short time-scale of the experiment. Our experimental design cap- 
tured a fraction of the potential effects of losing larger taxa on the 
functioning. We should expect that compensatory growth over the 
long-term would result in altered productivity and nutrient cycling. 
Further long-term investigations testing specifically for the conse- 
quences of changes in energy use with the loss of larger species would 
be necessary before the application of the current findings to broader 
realistic extinction scenarios. 

Global warming and overexploitation of resources are likely to 
induce a shift toward reduced species body size in aquatic ecosys- 
tems 18,19 , and will ultimately affect ecosystem functioning 20 . Studies 
linking directly body size to ecosystem functioning are essential to 
assess and predict subsequent effects of changes in size spectra on 
ecosystem dynamics. Our experimental study conducted in parallel 
in two aquatic systems showed that predictions of species loss effects, 
in multitrophic food webs, might not be efficiently achieved when 
based simply on body size, but showed great potential when com- 
bined with other easily measured ecological traits such as the trophic 
position. 

Methods 

Removal sequences. We ran two complementary experiments in parallel, with the 
same experimental design, both in marine and freshwater mesocosms 
(Supplementary Fig. SI). Starting from a regional pool of 9 marine and 10 freshwater 
taxa, we assembled the different mesocosm communities by removing 0, 1, 3, 6 and 0, 
1, 2, 4 taxa from the regional pool respectively in the marine and freshwater systems, 
according to two sequences: random and non-random. A non-random removal 
sequence denotes removals ordered by body size, from the largest to smallest. In the 
random sequences, all taxa were randomly removed from the regional pool such that 
each unique taxa was removed an equal number of times at each level in order to 
evenly distribute the effect of removing each taxa across the experiment (further 
details in ref. 17). All combinations were replicated three times, for a total of 93 
marine and 102 freshwater mesocosms. We did not perform any compensatory 
readjustment of biomass following taxa removal because the organisms we studied 
are discrete and have low population size, and it would have required too many 
assumptions about the redistribution of biomass among remaining groups (see 
discussion on this issue in ref. 17). The short duration of the experiment also 
prevented a population response to the removal of taxa. The experiment therefore has 
to be interpreted as a report of immediate effects of taxa removal on ecosystem 
functioning, before any compensatory growth takes place. 

Marine mesocosms. Marine mesocosms (21 L) were maintained during 6 weeks 
with a supply of filtered surface water ( — 50 um) from the Lower St. Lawrence 
estuary. Taxa consisted of 9 representative species of the neighbouring sublittoral: 
Cancer irwratus (70-80 mm), Strongylocentrotus droebachiensis (50-60 mm), 
Mytilus edulis (40-50 mm), Nucella lapillus (30-40 mm), Littorina littorea 
(20-30 mm), Gammarus spp. (10-20 mm), Testudinalia testudinalis (10-20 mm), 
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Figure 2 | Effects of taxa removal on ecosystem functioning according to random and non-random extinction sequences. The proxies of ecosystem 
functioning illustrated here are (a) periphyton dry mass (b) marine multifunctioning and (c) freshwater phytoplankton biomass. White boxplots indicate 
non-random extinction sequences and grey boxplots indicate random extinction sequences. The y-axis refer to the effect size between values of 
functioning recorded in mesocoms with removals and values of functioning recorded in reference mesocosms (no removal). 



Semibalanus balanoides (0-10 mm) and a mixture of Littorina saxatilis/obtusata 
(0-10 mm) (Supplementary Fig. S2). We chose length over body mass as potential 
measure of body size. Both approaches would induce biases by overestimating the 
metabolically active mass among organisms with heavy shells, compared with 
species of similar mass but having thinner shell, if any. In this experiment, length 
was preferred over body mass for technical issues, particularly in the freshwater 
system in which organisms exhibited very small body mass. Densities were adjusted 
to recorded densities where organisms were collected. We measured changes in 
periphyton dry mass, macroalgae and encrusted algae net production as proxies of 
ecosystem function. Periphyton dry mass was assessed from a 1 cm X 1 cm sample 
of pre-incubated Hester-dendy plate (dry mass; 24 hours at 60 J C). We estimated 
macroalgae net production by weighing the residual biomass from 100 g of the 
brown algae Fucus evanescens, placed in each of the mesocosms at the beginning of 
the experiment. We determined encrusted algae net production (Ralfsia verrucosa) 
in calculating changes in the algae cover present on a single rock placed in each 
mesocosm at the day 0, with image processing program ImageJ (National Institute 
of Health, USA). 

Freshwater mesocosms. Freshwater mesocosms (60 L) were filled with 40 L of 
filtered freshwater (20 um nylon mesh) and maintained outdoors for 8 weeks. Taxa 
consisted of small zooplankton (0.063-0.5 mm), large zooplankton >0.5 mm), 
Hyalellidae (3-4 mm), Dystiscidae (3-4 mm), Corixidae (5 mm), Coenagrionidae 
(6-14 mm), Planorbidae (6-11 mm), Gerridae (7-9 mm), Cyprinidae (22-32 mm) 
and Lymnaeidae (27-36 mm). We considered five ecosystem properties as proxies 
of ecosystem function: phytoplankton biomass, periphyton dry mass, 
bacterioplankton abundance, total nitrogen and total phosphorus concentration. 
Phytoplankton biomass was determined from 150 ml water samples filtered onto 
Whatman GF/F filters and extracted for 24 hours in 90% acetone, at 5°C in the 
dark 2 '. Chlorophyll a and phaeopigment concentrations were calculated after 
measuring fluorescence before and after acidification (HC1 1 M) with a 10-AU 
fluorometer (Turner Designs, Sunnyvale, USA) 22 . Periphyton dry mass was assessed 
as for the marine ecosystem. Bacterioplankton abundance was measured using 
standard flow cytometric analysis. Samples for bacteria abundance determination 
were fixed with glutaraldehyde 0.1% final concentration and stored at — 80 C until 
flow cytometry analysis 23 . Total nitrogen (TN) and total phosphorus (TP) were 
measured using the copper-cadmium standard reduction method for autoanalyzer 
after alkaline persulfate digestion 24 . 

Statistics. In all analyses, for the marine and the freshwater ecosystems, we 
considered respectively 3 and 5 measures of ecosystem functioning and a measure of 
multifunctioning. For multifunctioning, we conducted a PCA and used scores of the 
two axes as the dependent variable. The first and second axes explained respectively 
46.3% and 31.9% of the variation in the marine ecosystem (Supplementary Table Si). 
In the freshwater ecosystem, these axes explained respectively 25.3% and 22.5% of the 
variation. In the marine ecosystem, positive scores on the first axis were mostly linked 
to encrusted algae net production and negative scores to macroalgae net production. 
In the freshwater ecosystem, positive scores were attributable mostly to periphyton 
dry mass and negative scores to the combination of the four others functions 
measured. 

The effect of the body size of single taxon extinction on ecosystem functioning 
was tested with linear regression (Hypothesis 1). We also controlled with ANCOVA 
for the effect of the trophic position (grazers vs non-grazers) as a categorical 
variable. The difference between the effect of random and non-random extinction 
sequences was tested with an asymmetrical crossed two-way ANOVA, with 



Number of taxa removed and Sequence as main factors. The variability of the factor 
Sequence was then divided between: 1) random versus non-random extinction 
sequences and 2) among random extinctions sequences (Hypothesis 2, 
Supplementary Table S2). Normality was assessed using the Shapiro-Wilk test and 
homogeneity of variance was assessed using explanatory checks of plots of residuals 
against predicted values. 
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